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I. INTRODUCTION 

Our hope to cover the overpowering richness of Physics 
by microscopic theories constructed on simple elementary 
principles is based on the concept of universality, the pos- 
sibility of ignoring most of the short distance, microscopic 
parameters of the field theoretical models in describing 
the dynamics at finite scales [1]. Such an enormous sim- 
plification, obtained by inspecting the asymptotical UV 
scaling laws, is possible if the renormalized trajectory 
'spends' several orders of magnitude in the UV scaling 
regime and suppresses the sensitivity of the physics on 
the non-renormalizable, i.e. irrelevant parameters. 

But one tacitly assumes in this scenario that there 
are no other scaling laws in the theory or at least they 
do not influence the results obtained in the UV scaling 
regime. This is certainly a valid assumption for mod- 
els with explicit mass gap in the Lagrangian which ren- 
ders all non-Gaussian operator irrelevant at the tR fixed 
point. But what happens in theories with massless bare 
particles? Spontaneous symmetry breaking, or dynami- 
cal mass generations in general, may change the situation 
0, Q and open the possibility of tunable parameters of 
the theories which make their impact on the dynamics 
due to a competition between the UV and the IR scaling 
regimes. Such a competition can be highly non-trivial 
because both of these regimes can cover unlimited orders 
of magnitude of evolution for the RG flow and may offer 
a new way of looking at complex systems. 

The traditional way of 'understanding' a large system 
starts with the identification of its components and their 
elementary behaviour and then continues with the com- 
bination of these rules to come up with the whole picture. 
The competition between an UV and IR fixed points may 
introduce microscopical parameters which influence the 
dynamics in the macroscopic range only. Such a renor- 
malization microscope mechanism allows us to determine 
certain microscopical parameters by means of the long 
range dynamics and represents a microscopic, elemen- 
tary feature of a model which remains undetected unless 
the whole system is considered. 

We turn to a simple two-dimensional model in this 
work whose dynamics displays such a phenomenon. It 
has an explicit mass gap but this does not render the 



IR dynamics uninterestingly suppressed. This is be- 
cause the natural strength of a coupling constant is its 
value when expressed in units of the scale inherent of 
the phenomenon to be described, e.g. the gliding cut- 
off in the framework of the renormalization group (RG) 
in momentum space. Since the non-Gaussian coupling 
constants of the local potential have positive mass di- 
mension in two space-time dimensions, the usual slow- 
ing down of the evolution of dimensional coupling con- 
stants below the mass gap generates diverging dimen- 
sionless coupling constants, i.e. relevant parameters at 
the IR fixed point. But this is one part of the story 
only. The other part is that in order to have non-trivial 
competition between the UV and IR scaling regimes we 
need coupling constants which are non-renormalizable, 
i.e. irrelevant at the UV fixed point and turn into rel- 
evant at the IR end point. In order to render the cou- 
pling constants of the local potential non-renormalizable 
we need non-polynomial coupling. The simplest, treat- 
able case is the sine-Gordon (SG) model [4j whose La- 
grangian contains a periodic local potential and which 
possesses two phases. All coupling constants of the lo- 
cal periodic potential are non-renormalizable in one of 
the phases when the model is considered in the continu- 
ous space-time with non-periodic space-time derivatives 
We are thus lead to the massive sine-Gordon (MSG) 
model d, 0, [1, H, EH which has already been thoroughly 
investigated in the seventies as the bosonized version of 
the massive Schwinger model, I +1 dimensional QED, the 
simplest model possessing confining vacuum [|| . The dif- 
ferent UV and IR scaling laws have already been ad- 
dressed in this model @, H| but the careful identification 
of the full scaling regime, not only the asymptotical ends 
is needed to understand the global features of the RG flow 
diagram. This is the goal of the present work. One can 
only answer this question by solving the RG equations 
for the couplings and to find the effective potential. In 
this work we use the Wegner-Houghton (WH) functional 
RG method pd| , where the renormalization procedure is 
defined by the blocked action with gliding sharp cut-off 
to determine the flow of the Fourier amplitudes of the 
model. The handicap of the WH-RG method is that one 
cannot go beyond (if it even necessary) the local poten- 
tial approximation (LPA), since the gradient expansion 
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of the Lagrangian gives indeterminate evolution equation 
due to the sharp cut-off used. There are several other 
methods in the literature which are based on the evolu- 
tion of the effective action We show in this paper 
that the internal space RG method [13(, where the RG 
evolution is controlled by the mass and allows us in prin- 
ciple to go beyond the LPA gives the same results in LPA 
as the WH-RG method. 

Since the SG model has a condensed phase which arises 
at a finite cut-off fcsi one expects that if M < fcgi the con- 
densation also should appear in the MSG model. In our 
previous works [1, [3] we showed in the WH-RG frame- 
work that the non-trivial scalings appear in the deep IR 
limit for the SG and the MSG models, indeed. Since 
the SG model has a trivial, constant effective potential 
in either phase 114 . Hal . the RG methods based on the 
effective action [12|, LL3| may fail m treating both the SG 
and the MSG models. Nevertheless, when the control 
parameter approaches the to physical value of the mass, 
the internal-space RG enables one to find an evolution 
of the parameters which is analogous to their WH-RG 
flow. Also the sign of spinodal instability seems to ap- 
pear as a singularity in the internal-space evolution. In 
the latter case the analogy mentioned above allows one to 
change from the internal-space RG analysis to the WH- 
RG framework and to determine the phase structure of 
the MSG model. The situation seems to be similar as 
that for the WH-RG framework where the IR limit k — > 
of the blocked action is trivial but physics can be read 
off from its approaching this limit. 

A side-product of our RG analysis is that one recovers 
the well-known phase structure of the bosonized version 
of QED2 in the LPA. Therefore, one is lead to the con- 
clusion that the phase structure, as well as the spinodal 
instability should survive wave- function renormalization. 

Another possibility might be the usage of the Polchin- 
ski RG method which uses the blocked action and a 
smooth cut-off. Unlike the WH-RG method based on 
summing up the loop corrections during the evolution, 
Polchinski's procedure sums up the perturbation expan- 
sion, therefore it misses the spinodal instability in an 
obvious manner. In fact, one can easily show that the 
local potential evolved by this scheme is not of parabolic 
shape, a feature is thought to be a crucial sign of spinodal 
instabilities in LPA. 

The ionized (large (3 2 ) phase of the MSG model can 
be an example where thorough numerical investigations 
show that the bare UV irrelevant coupling becomes rele- 
vant in the deep IR regime. This fact clearly shows that 
the treatments of the MSG model based on any pertur- 
bation expansion is doomed to failure. 

The paper is organized as follows. In Sect.[TT]we derive 
the evolution of the coupling constants of the MSG model 
in the framework of the WH-RG method. In Sect. lIIII thc 
phase structure of the MSG model is presented. The 
RG microscope effect is discussed in Sect. lIVl The evolu- 
tion of the couplings constants in the internal-space RG 
framework is treated in Sect.lVIl We show in Sect.fVlthat 



our WH-RG results recover the well-known phase struc- 
ture of QED2, and finally, in Sect. lVIIl the conclusion is 
drawn up. 



II. BLOCKING IN MOMENTUM SPACE 

The MSG model is defined in 2-dimensional, infinite, 
Euclidean space-time by the Lagrangian 



S k = 



\{d^ x ) 2 + U k (cj )x ) 



(1) 



given in the leading order, local potential approximation 
(LPA) of the gradient expansion, k denotes the sharp 
momentum space cutoff and the potential is the sum 
Uk{4>) = \M 2 <j> 2 + Vk{4>), the second term being peri- 
odic, 



u n (k) cos(n(3 k (, 



(2) 



The blocking in momentum space [14|, the lowering of 
the cutoff, k — > k — Ak, consists of the splitting the field 
variable, <f> = (f> + <j>' in such a manner that <p and <p' 
contains Fourier modes with \p\ < k — Ak and k — Ak < 
\p\ < k, respectively and the integration over 0' leads to 
the WH equation |ll| 



(2 + kd k )U k {<f>) = - 



Ilnf 



(3) 



for the dimensionless local potential U k — k~ 2 U k . This 
equation is obtained by assuming the absence of instabil- 
ities for the modes around the cutoff. Only the WH-RG 
scheme which uses a sharp gliding cutoff can account for 
the spinodal instability, which appears when the restor- 
ing force acting on the field fluctuations to be eliminated 
is vanishing and the resulting condensate generates tree- 
level contributions to the evolution equation [10]. The 
saddle point for the blocking step, <fi' , is obtained by min- 
imizing the action, S k -Ak[<t>] = min0/ (S k [<j) + 0']) JEHU. 
The restriction of the minimisation for plane waves gives 



U k -& k {4>) = min 
p 



duU k {4> + 2pcos(iru)) 



. (4) 

in LPA where the minimum is sought for the amplitude 
p only. 

One can show that both evolution equations, Eqs. ([3]) 
and (|4]) , preserve the period length of the potential V k {<j>) 
and the non-periodic part of the potential, therefore 
M 2 = M 2 and (3 k — (3. Thus the mass is relevant pa- 
rameter of the LPA ansatz for all scales, 



Ml 



Ml 



(5) 



3 



A. Asymptotic scaling 

It is easy to find the asymptotic UV and IR scaling 
laws. One finds the evolution equation 

kd k u n {k) = ( - j 3 n -2 ) u n {k) (6) 

in the UV regime after ignoring 0(M 2 /k 2 ) and 
0{\U'^\ 2 /k 2 ) contributions with the solution 

MA) ( x ) { WT ^) ■ (7) 

The asymptotic IR scaling, well below the mass scale, 
is trivial because the mass gap freezes all scale depen- 
dence. The numerical solution of the complete evolution 
equation Eq. © is shown in Fig.Q] 

B. Impact of the mass gap 

It is instructive to compare the RG flow of the (mass- 
less) SG and the MSG models. The asymptotic UV evo- 
lution equations differ in 0(M 2 /k 2 ) terms only and the 
mass term gives small corrections to the scaling laws in 
this regime. But the mass gap freezes out the evolu- 
tion for any values of (3 thus more significant differences 
should show up between the SG and the MSG models. 
The dimensionful potential approaches a constant in the 
SG model as a result of the loop-generated or instability 
driven evolution in the ionized or the molecular phase, re- 
spectively [3, [H| . The evolution of the potential freezes 
out below the mass gap of the MSG model and a non- 
trivial potential is left over in the IR end point, reflecting 
the state of affairs at k w M. The IR scaling is trivial for 
k < M, u n (k) ~ k~ 2 . In order to go beyond the asymp- 
totic scaling analysis and to find out more precisely the 
changes brought by the non-periodical mass term to the 
RG flow we have to rely on the numerical solutions of the 
evolution equations. 

Let us consider first the regime f3 2 > 8ir which is free of 
spinodal instabilities in the massless case. The evolution 
of the first four coupling constants, u\, . . . , U4 is shown 
in Fig. m for 1 > 8ir. The UV scaling regime is confined 
in this plot to the very beginning, around k/A f=s 1 [Til ], 
and what we see here is that the flows of the SG and 
the MSG models agree for 1 > 8ir even in the IR scal- 
ing regime down to the mass gap. The freeze-out below 
the mass gap takes place naturally without generating 
spinodal instabilities. 

The comparison of the massive and the massless cases 
is more involved for (3 2 < 8ir due to the appearance of 
instabilities. If the scale fcgi where instabilities appear is 
higher than the mass gap, ksi > M, then the RG flows 
of the MSG and SG models are similar down to M and 
they, i.e. both display instabilities and differ for k < M 
only, as shown in Fig.O When the freeze-out scale is 
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FIG. 1: RG flow of iti, . . . ,u 4 for p 2 = 12tt and M 2 = 1(T 9 A 2 
(solid line) or M 2 = (dense dotted line). The plot also 
shows the asymptotic IR scaling ~ k~ 2 (rare dotted line). 
The massive and the massless flows depart at the scale M, 
indicated by a dashed vertical line, placed at the intersection 
of the dotted lines. 
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FIG. 2: RG flow of ui, . . . ,M 4 for f3 2 = 4tt and M 2 = 10~ 9 A 2 
(solid line) or M 2 = (dotted line). The scales of the spinodal 
instability fcsi/A and the mass M/A are indicated. The inset 
shows RG flows of u\ for several different initial values at 
M 2 = 10" 9 A 2 . 



reached first during the evolution, i.e. the scale ksi of 
the SG model with the same potential as V(cj>) of the 
MSG model is smaller than M then the instability does 
not occur in the MSG model. 

It is instructive to determine the boundary of the re- 
gion in the coupling constant space with spinodal insta- 
bility with the truncation where a single Fourier mode is 
kept only, u n = u5i n . The condensate appears during 
the evolution at scale fcsi, satisfying k^+M 2 + U{{ (<f>) = 
for some <fr. The approximate analytic expression, 
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Eq. 0, gives 



10" 



(A 2 + M 2 ) 



A 2 + M 2 
/3 2 w(A) 



Af 2 



(8) 



for /3 2 < 87r, suggesting that the coupling constant can 
be weak enough to allow the mass term to remove the 
condensate. The RG flow obtained numerically for u±(k) 
shown in the inset of Fig. [2] confirms, as well, that the 
mass can be strong enough to prevent the formation of in- 
stabilities. One can get an estimate of the critical value of 
the coupling constant by equating kgi and M in Eq. ([8]), 



uic(A) = 



A 2 + M 2 

/3 2 



2A/ 2 



A 2 + M 2 



(9) 



The boundary of the region with instability is shown in 
Fig. [3] In contrast to the SG model where the instability 
extends over the whole phase with 1 < 8ir the mass 
term always wins at the IR end point of the flow of the 
MSG model and removes the condensate at some low but 
finite value of the scale k. 

The disappearance of spinodal instability and the triv- 
ial scaling, |{t„| ~ k~ 2 7] can be made plausible also 
by the following, simple analytic consideration. Namely, 
Eq. ([3]) can be rewritten as 



^B niS kd k u s (k) = ^2 

s s 

with 
resulting 



a-ifi n k 
k 2 + M 2 



u 8 (k), 
(10) 



s/3 2 (n-s) 2 U|„_ s |(0)-(n+.s) 2 w n+s (0) 



2(k 2 + M 2 ) 



s(3 2 a n 



2M 2 



(11) 



(12) 



an RG invariant quantity in the IR scaling region. Fi- 
nally, the second term on the right hand side of Eq. (flTJ|) 
can be neglected for k <C M, yielding u n (k) ~ u n (0)k~ 2 
in the IR scaling region. 



III. PHASE STRUCTURE 

The mass term deforms the phase boundary of the 
SG model by extending the ionized phase. In this 
phase of the SG model the IR scaling law generates 
the scale dependence of the coupling constants u n (k) 
with n > 2 through u\(k), namely renders the ratios 
i?f G = u n (k)/ui(k) RG invariant. It was checked nu- 
merically that Rn ISG = \u n (k)\ / Ui (k) is RG invariant in 
the IR scaling region of the MSG model without conden- 
sate and the potential at k = depends on the initial 
value of ui only. 
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FIG. 3: The region with spinodal instability. Different sym- 
bols show the boundary obtained by the numerical computa- 
tion with 10 coupling constants and the lines show the results 
obtained by taking into account a single Fourier mode only. 



The phase with condensate is similar to those of the 
SG model. The potential develops quickly into a super- 
universal, initial condition independent shape [Tij when 
M < k « fcsi, cf. the inset of Fig.fSJ But this scaling 
regime ends at k ss M where trivial scaling laws come 
into force down to k = 0. The matching of the IR scaling 
of the SG model [l4| with the trivial scaling law gives 
Un (0) = (~l) n+1 2M 2 /n 2 (3 2 . 

The modification of the phase boundary induced by 
the mass can be seen by means of the sensitivity matrix 
[H,[l4|], too. This matrix, defined as the derivatives of the 
running coupling constants with respect to the bare one, 



du„(k) 



(13) 



™' m du m (A) ' 
develops singularities when the UV and IR cutoffs are 
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removed at the phase boundaries only. The typical be- 
haviour is depicted in Fig.dJ showing that the appear- 
ance of the condensate generates first singular turns and 
leads later to radically different scale-dependence in this 
matrix. 
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FIG. 4: The (1,1) element sensitivity matrix is computed 
numerically for various initial conditions for /3 2 =47T as the 
function of k/A. The trajectories of molecular and the ionized 
phases are shown by dashed and solid lines, respectively. 

The molecular phase of the two-dimensional SG model 
was compared with the four-dimensional Yang-Mills the- 
ory [l4| because both support spinodal instabilities in 
their vacuum and the periodic symmetry of the SG model 
is formally similar to the center symmetry of the gauge 
theory. This analogy can be extended to the MSG model, 
where the breakdown of the periodicity which manifests 
itself in the scaling laws for k < M is comparable with the 
spontaneous symmetry breaking in a four dimensional 
Yang-Mills-Higgs system. The solitons cease to be sta- 
ble and elementary plane- wave excitations appear in the 
scattering matrix of the MSG model due to the break- 
down of the periodicity by the mass term. In a similar 
manner, the confining forces are lost by the spontaneous 
breakdown of the global symmetry and elementary ex- 
citations propagate in the vacuum of gauge theories. In 
other words, the competition of the confinement and the 
symmetry breaking scales in the gauge theory is similar 
to what happens between kgi and M in the MSG model. 



IV. RG MICROSCOPE 

The non-triviality of the IR scaling regime opens the 
possibility of having relevant non-renormalizable oper- 
ators, being irrelevant around the UV fixed point, but 
becoming relevant in the IR scaling regime. To find out 
whether this can be the case let us take the running cou- 



pling constants, given by Eq. |(7J 

u n (k) = (-l)" +1 <(A)i?^ SG fc 2n - 2 

fk\-' 2n fk 2 + M 2 \^ , . 

X U) Ut^J ' (14) 

and calculate the (1,1) element of the sensitivity matrix 
in Eq. (flU)) in the limits k — > and A — > oo, 

S m ~AT 2 At. (15) 

We set /3 2 > 8ir, i.e. move deep in the ionized phase in the 
language of the SG model where all coupling constants 
are non-renormalizable. But the IR scaling law makes 
the coupling constants relevant with its factor k~ 2 . More 
precisely, the dimension 2 coupling constants freeze out 
meaning that their values counted in the 'natural' units 
of the running cutoff, k, diverge as the IR end point, 
k = 0, is approached. The result is the regain of the sen- 
sitivity of the long distance physics on the choice of the 
bare, microscopic, non-renormalizable coupling constants 
which was suppressed in the UV scaling regime. As long 
as the UV cutoff A and the IR observational scale k are 
chosen according to k = t~ Pln /i and A = t Pvv fi where 
Pir = .Puv(/3 2 /87r — 1) > (with any t > 1 and fx > 0), 
the observed IR dynamics depends on the choice of the 
bare coupling constant. This rearrangement is reminis- 
cent of the traditional microscope in the sense that the 
amplification, i.e. divergence of the RG flow in the IR 
scaling regime balances the suppression, i.e. focusing of 
the flow in the UV scaling regime. As a result the large 
scale observations can fix the value of the microscopic 
parameter with good accuracy. Notice that this kind of 
dynamics must be put in the theory at microscopic scale 
even though it starts to influence the physics in the IR 
scaling regime. Such an interplay between the different 
scaling regimes represents a highly non-trivial, global ex- 
tension of the simple universality idea which is based on 
the local analysis of the RG flow at a given fixed point 
and might be a key to phenomena like superconductivity 
and quark confinement 0, Q ■ 



V. TWO PHASES OF QED 2 

The most remarkable consequence of the non-trivial 
phase structure of the MSG model is that the two- 
dimensional quantum electrodynamics (QED2) also has 
two phases. The Lagrangian of the QED2 is given as 

C = -\f^F^ + ~ ieA M )^ " (16) 

where — d^A v — d v A^, m and e are the bare rest 
mass of the electron and the bare coupling constant, re- 
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spectively. The bosonization rules gives 



: iptp 
in =■ "07m 



-cmM cos(2v / 7T0), 
-cmM sin(2v / 7r0), 

-=eu V d v 4>, 

A 

-N m (d^)\ 



(17) 



where N m denotes normal ordering with respect to the 
fermion mass m, c= exp ( < y)/2n with the Euler constant 
7, and M = e/y/"K the 'meson' mass. The MSG Hamil- 
tonian 



H = N 



M 



1 



-it 



i 



{di<j) x Y + -M Z 4> A X + ui cos(/3<fc) 



for /3 2 = 4-7T and iti = cmM is the bosonized version of 
the QED2 . According to Eq. (O the dimensionful electric 
charge e = M^pH does not evolve. Using Eq. ([7]) the flow 
of the electron mass m(k) 



m(k) — m(A) 



fk 2 



\A 2 



M 2 



m(A) 



mi(A)v^F 



(19) 

inherits all the properties of the flow of the first Fourier 
amplitude u\(k). It implies that there exists a critical 
value 



m c (A) 



3 (A)0F V2 



47TC 



(A 2 + e 2 /^) 1 / 2 



(20) 



of the bare mass which determines whether the IR value 
of the mass depends on its UV value or not. For m(A) > 
m c (A) the evolution runs into the spinodal instability 
and the IR value of m becomes independent of its bare 
value: 



m(k -> 0) 



2 ev ^F 
c4tt 2 ' 



(21) 



while for m(A) < m c (A) Eq. ([19)) implies the IR be- 
haviour 



m(k -> 0) = m(A) 



M 2 



A 2 + M 2 



1/2 



(22) 



Therefore the sensitivity matrix in the parameter space 
(to, e) taking the values S± t i = for m(A) > m c (A) and 



= 



dm(k) 



k 2 + e 2 /ir 



dm(A) \A 2 + e 2 /n 



1/2 



(23) 



for m(A) < m c (A) indicates the existence of two differ- 
ent phases in QED2. Lattice calculations also affirmed 
this result [l?], [3- For large coupling (e >> to), the 
model has a unique vacuum at ip = 0. For weak coupling 
(e << to), the reflection symmetry is spontaneously bro- 
ken and the model has non-trivial vacua, located approx- 
imately at (j> = ±y / 7t/2. According to lattice simulations 



and density matrix RG studies of the MSG model the 
critical value which separates the two phases of the model 
is TO/e c = 0.3335. The analytical result m/e c — 0.3168 
for the critical point in Eq. (|2U)) suggests that the RG 
methods using LPA enables us to determine the phase 
structure of the MSG model in a reliable manner. 



VI. INTERNAL SPACE RENORMALIZATION 

One can go beyond the LPA by using the internal space 
RG method. We define the generating functional of the 
connected Green functions as 



W\j] = log J V<h 



-s B m+3-<i> 



(24) 



with external source j x . The shorthand notation / • g = 
J x fx9x is used. The effective action is defined as the 
Legendre-transform of W[j], 



r[<p] = j><p-w\j\, 



(25) 



where the external source j x is expressed in terms of Lp x 
according to the implicit equation 



SW[j] 
Sj ' 



(26) 



The idea of internal-space RG is to eliminate quantum 
fluctuations successively ordering them according to their 
amplitudes. This can be achieved by introducing an ad- 
ditional mass term into the action, 



SxW = S B [ 



1 



2 J.2 



-X 4> 



(27) 



with the control parameter A. For A = Ao being of the 
order of the UV cut-off A the large-amplitude fluctuations 
are suppressed and decreasing the evolution parameter A 
towards zero, they are continuously accounted for. Let us 
separate off the suppressing mass term from the evolving 
effective action, 



r x [<p] = f A M + lAV, 



and use the ansatz 



±(d^ x ) 2 + u x (<p x ) 



(28) 



(29) 



with U x (ip) = \M\lp 2 + V x {<p) and V x (<p) = 
YlnLi u n(X) cos(n(3\ip). The functional evolution equa- 
tion is 



1. 



-Tr 



\2 r , p(2) 
°x,V + 1 A x,y 



where T 



(2) 

A x,y 



6 2 T x /6(p x 5(p y , ll|, reads 



p 2 + A 2 + M 2 + VI 1 (ip) 



(30) 



(31) 



7 



for homogeneous field configurations ip for the po- 
tential V\(cp) of the ansatz Eq. ([25)1 where V x (ip) = 
d 2 V\(ip)/dip 2 . The two-dimensional momentum integral 
can easily be performed, giving 



(i+x 2 d X 2)v x = — log 

07T 



(A/A) 2 + 1 + Ml + 



VI' 



Ml 



(32) 



where the dependences on the field variable tp is sup- 
pressed. By performing the Fourier expansion in both 
sides of Eq. (|32|) one obtains a set of differential equa- 
tions for the couplings « n (A), /3 X and AI X . Since the left 
hand side of Eq. (|3ip does not contain polynomial terms, 
the mass parameter does not evolve, M 2 = M 2 . Thus 
the mass is a relevant parameter of the LPA ansatz for 
all scales, 



M 2 X =\- 



l M 2 



(33) 



cf. Eq. (O. Furthermore, as in the case of the WH-RG 
equations, if the right hand side at a given scale A has 
a given period 2ir/[3 x in the internal space, then V\-d\ 
shall have the same period, so that the parameter (3 X = (3 
evolves neither. Note that the argument of the logarithm 
on the r.h.s. of Eq. (f3"2")) should be positive, otherwise the 
evolution fails. This criterion is reminiscent of the WH- 
RG method and it may generate spinodal instability like 
singularities on the internal-space RG renormalized tra- 
jectory, to be interpreted as quantum phase transition. 
In order to obtain the evolution equations for the cou- 
plings more easily we differentiate Eq. ([32]) with respect 
to tp, 

1 A 2 ~ 
- ——VI" = 
8tt A 2 a 

A 2 
A 2 



+ 1 + M 2 + Vl' ) (1 + M 2 + Vl')(l + X 2 d X 2)Vi. 

(34) 



e>(A 2 ) > m 2 > 



/ A I . At these 'UV scales the large-amplitude field fluc- 



The evolution should be started at Aq 

tuations with |<y9| 2 ^> |V^'|/A 2 are suppressed and the 
effective action can be calculated perturbatively. 

A. Asymptotic scaling 

The asymptotic scaling for A 2 satisfying A 2 ^> A 2 ^> 
M 2 3> \Vl'\ can be deduced by using the independent 
mode approximation, 



f 



(l + \'d X 2)V x = - — 
8n 



1 



1 + M 2 A 2 /\ 2 + l + M 2 



V", 

where the field independent terms are neglected. The 
evolution of the couping constants decouple and one finds 



"„(A! = r/„(A ll ) ( j 



A 2 + M 2 



A 2 + A 2 + M 2 



(36) 



for the dimensionless couplings constants u n (A), yielding 
u n ~ A" 2 ^ 2 / 4 ^ 2 . (37) 

B. Ionized phase 

In order to find out the non-asymptotical regime one 
has to solve a system of coupled evolution equations for 
the couplings constants u n (\) as the control parameter A 
is decreased from Ao = A down to A = 0. The evolution 
equations has been derived and solved for (3 2 > 8tt nu- 
merically. It was found that the increase of the number 
of the couplings u n beyond n = 10 does not influence 
the evolution of the first few couplings, similarly to the 
WH-RG equations. 

The evolution of the first four coupling constants, 
£ti,...,U4 is shown in Fig. [5] for /3 2 = 12-7T. The inter- 
nal space RG method gives qualitatively the same scaling 
laws both in the 'UV and in the TR' regions just as the 
WH-RG method. The numerical value of the fundamen- 
tal coupling constant ui(A) was found to follow closely 
the analytic form of Eq. (|36[) . The renormalized trajec- 
tory shares the feature, known from the WH-RG scheme 
that the SG and the MSG models with 1 > 8ir agree in 
the IR scaling regime down to the mass gap. 



1. IR scaling 

The decrease of the control parameter A drives us out 
from the asymptotic region. According to Fig.[S] the IR 
scaling is u n ~ A™^ /4tt-2) , Such a scaling behaviour can 
be obtained analytically from the functional RG equation 
in Eq. (f3"2")l which becomes 

(1 + \ 2 d X 2)Vi + V{'(1 + \ 2 d X 2)Vi = -i-Vf (38) 



for M 2 < A 2 < A 2 . Making the ansatz 
u n {\)=c n \ nv 



(39) 



with r) > 0. For n — 1 one gets r\ = /3 2 /47r — 2 > 0. For 
n > 1 we obtain the recursion relation 



\P 2 S"=i(2 + sy)s(n - s) 2 c n . s c s 
n(2 + nrj- n 2 §^) 



(40) 



The coefficients c„ can be expressed in terms of c\ , since 
ci = ui(A)(A/A)'> and therefore c„ = (-l)" +1 ^(A)i? ri , 
with R\ — 1 and all R n being independent of the bare 
couplings. These properties were confirmed numerically. 
They imply that the dimensionless effective action in the 
IR regime can be parametrized by the single bare param- 
eter Wi(A). 

In order to consider the effect of the mass M in the 
theory we also determined the RG evolution of the cou- 
plings for the (massless) SG model, which depicted by 
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FIG. 5: RG flow of Mi, . . . ,u 4 for f3 2 = 12vr and M 2 = 1(T 9 A 2 
(solid line) or M 2 = (dense dotted line). The plot also 
shows the asymptotic IR scaling ~ A~ 2 (rare dotted line). 
The massive and the massless flows depart at the scale M, 
indicated by a dashed vertical line, placed at the intersection 
of the dotted lines. 



dense dotted lines in Fig.[5l One finds that the scaling 
behaviour in Eq. (|39[) goes on to infinitesimal values of A. 
We note that all couplings are irrelevant and the effective 
action is zero in the A — > limit. After introducing the 
mass M it is clear from Fig. [5] that the evolution of the 
potential freezes out below the mass gap and a non-trivial 
dimensionless potential is left over in the IR end point. 
The IR scaling is trivial for A < M, u n (X) ~ A~ 2 . Since 
the evolutions for the SG and MSG models are identical 
down to the scale A ~ M, the evolving effective action of 
the MSG model inherits the properties of the SG model, 
namely it depends on the initial value of the fundamental 
mode ui(Xo) only. In fact, it was found numerically that 
Rn = \ u nW\/ u iW 1S KG invariant in the scaling 
region A < M. 

One sees also that for 1 > 8tt theories with various 
values of the mass M belong to the same phase. This 
is conclusion arises because one detects the same scal- 
ing behaviour of the dimensionless parameters u n (X) for 
X/M < 1 and the same qualitative behaviour of the RG 
invariant constants R¥ SG for all values of M. 



C. Molecular phase 

The typical evolution is depicted for several bare val- 
ues iii (A) in Fig. [5] for 1 < 8tt. Numerics shows im- 
mediately, that the evolution stops at a non-vanishing 
value of A = A c due to the appearing a negative argu- 
ment of the logarithm in Eq. ([32]) for the Ux(A) > u\ Cl 
where ii\ c is a critical value of the first Fourier amplitude 
of the bare periodic potential. Then the evolution equa- 
tion loses its validity and presumably an alternative RG 



FIG. 6: Internal-space RG evolution of the coupling u\ for 



several different initial values at M 



l(T a A a 



equation is necessary as in the WH-RG framework, where 
the appearance of the spinodal instability implies a tree- 
level blocking relation [H, 0, [l5| . Though plausible but 
is not obvious that the singularity in the internal-space 
evolution is also rooted in the spinodal instability. The 
clarification of this point needs further efforts. 

For bare values iii(A) < u\ c the evolution does not 
stop and goes below the mass scale. Then, as in the case 
of (3 2 > 8ir, the trivial scaling u n (X) ~ A~ 2 is obtained 
in the limit A — > 0. 

Although one cannot follow the RG evolution for A < 
A c for iii (A) > U\ c by means of Eq. (|32p . one can still 
determine the critical value u\ c analytically, when the 
local periodic potential is restricted to its first Fourier 
mode. The singularity appears during the evolution at 
the scale A c satisfying X 2 -, + M 2 + V£ (<p) = for some 
if. Using Eq. ([55)) as a good approximation of the scaling 
for A > A r , one finds 



A 2 = A 2 (/?V(A)) 



— M 



(41) 



for (3 2 < 8tt. The negativity of the right hand side sug- 
gests that the coupling constant is sufficiently weak to 
allow the mass term to remove the singularity. For the 
opposite case X 2 > one can estimate the critical value 
of the coupling constant by equating A c to M in Eq. (I41|) , 



1 /2ikP 



A 2 



(42) 



The value of A c is plotted on the plane (f3 2 /tt,ui(A)) in 
Fig. [71 In contrast to the SG model where the singularity 
appears for all i*i(A) for 1 < 8tt the mass term always 
wins at the IR end point of the evolution for the MSG 
model and removes the singularity at some low but finite 
scale A r . 



9 




10 u 
10" 1 

10~ 3 

10- 4 



FIG. 7: The scale A c , where the singularity appears, as the 
function of the initial value of the first Fourier amplitude ui c 



and f3y-K for M 



10" 



The introduction of the wave function renormalization 
might modify the phase structure. Nevertheless, one ex- 
pects that both 'molecular' phases survive wave-function 
renormalization. This is because the MSG model with 
the particular value of 1 = Air represents the bosonized 
version of the 2-dimensional quantum electrodynamics 
QED 2 which exhibits two phases. 



the model possesses condensate of elementary excitations 
with non-vanishing momentum, spinodal instability, for 
weak enough mass in the remnant of the molecular phase 
of the SG model. This condensate, the sign of the period- 
icity of the local potential, generates non-trivial effective 
potential and phase structure despite the explicit, stable 
mass for elementary excitations in the deep IR region. 
The sensitivity matrix allows us to study the way the 
ultraviolet parameters influence the IR physics. It was 
found that the suppression of the sensitivity on the non- 
renormalizable bare coupling constants, generated in the 
UV scaling regime, can be overturned by the increasing 
sensitivity piled up in the IR scaling regime if the UV 
and the IR cutoffs are removed in a coordinated manner. 
As a result, a non-trivial, global extension of the univer- 
sality is found which goes beyond the local studies of the 
RG flow around the UV fixed point only. 
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